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Abstract. In a clean Fermi liquid, due to spin up/spin down symmetry, the dc spin 
current driven by a magnetic field gradient is finite even in the absence of impurities. 

Hence, the spin conductivity (7g assumes a well-defined collision-dominated value in 
the disorder-free limit, providing a direct measure for the inverse strength of electron- 
electron interactions. In neutral graphene, with Fermi energy at the Dirac point, the 
Coulomb interactions remain unusually strong, such that the inelastic scattering rate 
comes close to a conjectured upper boimd t^^i ^ kBT/h, similarly as in strongly 
coupled quantum critical systems. The strong scattering is reflected by a minimum 
of the spin conductivity at the Dirac point, where it reaches ag = ^^"^jt weak 
Coulomb coupling a, ~ A*-B being the magnetic moment of the electronic spins. Up 
to the replacement of quantum units, e^/fi — > /ig/fi, this result equals the collision- 
dominated electrical conductivity obtained previously. This accidental symmetry is, 
however, broken to higher orders in the interaction strength. For gated graphene, and 
2d metals in general, we show that the transport time is parametrically smaller than 
the collision time. We exploit this to compute the collision-limited analytically as 
as = ^ (^)^ ^ with C — A-'K^a^ [| ln(l/2a) — l] for weak Coulomb coupling a. 



PACS numbers: 72.25.Rb, 67.90.-hz,71.10.-w,73.23.-b 
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1. Graphene - a strongly coupled, relativistic electron-hole plasma 

Graphene is a monolayer of graphite, which forms a zero-gap semiconductor whose low 
energy quasip articles obey the massless Dirac equation [H El [3] . At charge neutrality, 
the Fermi surface reduces to two inequivalent Fermi points, forming a non-analyticity in 
the density of states, which can be viewed as a very simple quantum critical point jl]. 
However, on top of that, as a consequence of the vanishing density of states and the 
linear dispersion of the 2d quasiparticles Coulomb interactions remain unusually strong. 
They are only marginally irrelevant under renormalization, flowing only logarithmically 
to zero with decreasing temperature T [5]. Moreover, screening is suppressed so that 
the long range character of the Coulomb interactions remains mostly intact, except for 
screening due to thermally excited carriers. This strong interaction is reflected, e.g., 
in the inelastic transport scattering rate being proportional to Na'^T, where = 4 is 
the number of fermionic species, and we use units with ks = 1. Here, a = e^/nhvp 
is the dimensionless "fine structure constant" characterizing the strength of Coulomb 
interactions, where k is the average dielectric constant of the adjacent medium and 
fi? is the Fermi velocity of the linearly dispersing quasiparticles. The large scattering 
rate in graphene comes close to saturating a kind of Heisenberg uncertainty principle 
for quasiparticles |6]. This latter is the conjecture that the scattering rate can never 
significantly exceed the thermal energy scale, hr^-y < T. Upon approaching this limit by 
dialing up a — )■ 0(1), one indeed expects to drive a quantum phase transition towards 
an insulator [3 |H1 with completely different low energy excitations. 

Due to the strong marginal interactions, the neutrality point of graphene exhibits 
a transport phenomenology very similar to quantum critical points in more complex, 
strongly coupled materials [U HDl [HI [12], [13]. What makes graphene particularly 
attractive in this context is the fact that it is probably the simplest condensed matter 
system to possess the ingredients for this strong coupling phenomenology, being at the 
same time rather easy to produce experimentally at a high level of purity [15] . 
Indeed, a spectacular experimental proof for the non-negligible Coulomb interactions in 
clean graphene is the recent observation of the fractional quantum Hall effect [Ml |T7l HE] . 

The eventual logarithmic flow of the Coulomb coupling constant a toward zero 
justifies a perturbative analysis of the transport properties at low enough temperatures. 
Upon extrapolation to moderate coupling this captures the gist of the transport 
phenomenology, even though the results should not be trusted at a quantitative level 
beyond the weak coupling regime. In the "quantum critical window", i.e., at small 
chemical potential of the carriers, \^\ < T, the quasiparticles of graphene form an 
interacting "hot" (non-degenerate) electron-hole plasma with rather unusual transport 
properties, such as an anomalously low viscosity and a concomitant tendency toward 
turbulent electronic current flow in mesoscopic samples. [19] Those are features one 
should equally well expect in strongly coupled quantum critical matter, if it is sufficiently 
clean. Another interesting consequence of the strong scattering in graphene is the 
relativistic hydrodynamic transport which emerges in a regime of moderately low 
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frequencies below the inelastic collision rate. [201 |2T] 

At finite carrier density, a simple estimate of the inelastic scattering rate in random 
phase and Born approximation suggests that 

^^ne\-max(T,H)^^^-^^, (1) 

where a ~ a(e) = aA/[l + a;A/41n(A/e)] denotes the renormalized strength of Coulomb 
interactions a at the relevant energy scale e = max[/i,T], whereby A is a UV cutoff [5]. 
(In the sequel we simply write a with the understanding that this renormalized value 
should be used.) At finite /x, the scattering rate decreases rather quickly with T, 
following the familiar law T^/|/i|, which we will show to be merely logarithmically 
dependent on the interaction strength a. 

As was first pointed out in the context of the superfluid-insulator quantum phase 
transition [12] the particle-hole symmetric point /i = exhibits a finite collision- 
dominated conductivity, even in the absence of impurities. Indeed, the application 
of an external electrical field induces counter propagating particle and hole currents, 
and thus no net momentum. The latter is usually the source of infinite current response 
unless the momentum decays due to impurities, Umklapp scattering being negligible 
at the temperatures we have in mind. However, in neutral graphene one finds a 
disorder-independent conductivity which is solely due to electron-hole friction. Scaling 
arguments [221 IIH] based on the Drude formula using the thermal density of carriers 
nth ~ {T/hvpY, the above discussed inelastic scattering rate, and a T-dependent 
effective mass rrics ~ T/vp suggest a conductivity which is nearly constant, apart from 
a weak logarithmic growth due to the renormalization of a, 

2 2 

ae/i = = — — — . (2) 

nies a{T)'^ n 

This is indeed confirmed by a microscopic calculation based on the semiclassical 
Boltzmann equation, which becomes asymptotically exact for T ^ A where a -C 1, 
yielding the pref actor k = 0.121 [in], [f] 

A similar phenomenon arises when we consider spin transport in the presence of 
magnetic field gradients. A finite, collision-limited spin conductivity cr^ or spin diffusion 
constant Ds arises whenever the symmetry between spin up and spin down is not broken 
by a background field. Unlike the electrical conductivity, which diverges at finite gating 
in the disorder free-case, spin diffusion may serve as an interesting and direct measure 
of electron-electron interactions in clean metals in general. In this work we analyze 
spin transport in the linear response regime as a function of gating. This is motivated 

I This weak coupling approach could be extended to include screening effects along the lines of the 
recent work [23", by considering the simultaneous limit a — > and N ^ oo, keeping aN constant. 
In this limit the random phase approximation becomes exact and could be incorporated in the kinetic 
equation. The above calculation, as well as the one in the present paper, correspond to the limit aN <^ 1 
of such an approach. On the other hand, for aN ^ 1 one will find a scattering rate fi-Tj^^i ^ T/N and 
a conductivity Ue ~ while remaining in a regime where the perturbative weak coupling approach 
is still justified. 
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Figure 1. Finite collision-dominated spin transport. In the presence of a magnetic 
field gradient spin down carriers are accelerated opposite to spin up carriers, creating 
a counterflow with vanishing total momentum. Interactions among the carriers lead to 
friction between the counter propagating up and down currents, which degrades the 
spin current and leads to a purely interaction-limited spin conductivity. 



in particular by recent measurements of the gating dependence of spin diffusion in 
graphene [2H [2S], which found the spin diffusion constant to be roughly proportional 
to the charge diffusion constant. In the relatively disordered samples under study, this 
proportionality was attributed to impurities which presumably affect the mean free 
path of both diffusion processes in a similar way. However, this proportionality is rather 
unlikely to survive in clean, suspended graphene where interactions are expected to 
become the dominant source of scattering. Clearly the simple proportionality cannot 
hold at finite gating where o"e remains sensitive to impurities, while as is entirely collision 
dominated in the weak disorder limit. As we will see, also in the ungated case, at the 
Dirac point, the comparison of the two diffusion constants may provide an interesting 
probe for strong coupling effects. 



2. Model 



At low energies, graphene is described by massless Dirac quasiparticles obeying [3] 



. , ^ L ^ . . _ I L ^ 



! ^ 'ucrXk 'i^o-Xk^ 



(3) 



where vp is the Fermi velocity and are the Pauli matrices in sublattice space. The 
indices {z/, a, /}, label the valley {u = ± for K or K' point, respectively), the electron 
spin, and the valence or conduction sub-band, respectively. The quasiparticle operators 
^i/o-ifc diagonalize the Hamiltonian. 

In this basis, electron-electron interactions are described by [10] 

d^fci d^fca d^g rp (h h ^ ^ t 

(27r)2 (27r)2 (27r)2 ^'^'^'^'^^' ^'^^^*A4fc4^iA3fc3^i^2fe^^^i 



Hi 



ij A1A2A3A4 

where z,j stand for index pairs (z^, cr), and 

V{q) 



T. 



A1A2A3A. 



1 + A4A] 



l^4| li^ll 



, (5) 
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V{q) being the Fourier transform of the two particle interaction. Here and in the 
following we will use the notations = — q and ^4 = fci + q, as well as Ki = kf + ikf. 
Below we will mostly consider Coulomb interactions, 

^Cb(?) = — = (6) 

where the "fine structure constant" a has a slow logarithmic flow under RG 
transformation as mentioned above. We will be mostly concerned with a weak coupling 
analysis to leading order in a. Close to the Dirac point, fi < T, the inclusion of 
screening yields only subdominant corrections and has thus not been included in earlier 
studies on collision-limited transport in that regime flUl [26] . However, in the gated 
case, fi ^ T, the screening of long range interactions limits the scattering rate at small 
momentum transfers. Since our results do not depend crucially on the form of this 
cut-off, we content ourselves with a static (Thomas-Fermi) screening approximation of 
the Coulomb interaction 

V{q) = = (7) 

Here go = SaT ln[2 cosh(/i/2T)] is the Thomas-Fermi wave vector, which tends to 
go = 4a/x in the limit /i ^ T. A detailed analysis of screening effects at the Dirac 
point, including dynamical effects, has been given in the recent work [23]. [|] 

2.1. Spin current 



Here we are interested in the spin current driven by a magnetic field gradient. The dc 
magnetization current density due to moving electronic spins is given by 



where fig ~ /^s is the magnetic moment of the electrons, = X^vp is the group velocity, 
and f^^ij^ is the occupation number of the electronic quasiparticles. For a discussion of 
coherent contributions to the current at finite frequency, see Ref. [TT]. In linear response 
to a gradient in an external Zeeman field V-B (perpendicular to the graphene plane, see 
Fig. [1]) we write for the latter 

= fxk + g^xkfxkil - fxk)f^s\VBl (9) 
where f^f^ = l/(e[^^^'^'^~'^)/^ -|- 1] is the equilibrium distribution function, and 

9aXk = ^ik-^9ax{k), with e = ^— . (10) 

I Vi? I 

§ We note that the approximate analysis of coUision-hmited transport in Ref. [23^ is in quaUtative 
agreement with previous calculations |10( I20j and our results here, which, however, yield the exact 
leading order in a. 
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Since / and g do not depend on the valley index u, it is suppressed here and below 
Finally, the spin conductivity is defined as 



\VB\ 



2E 



(2vr) 
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2.2. Boltzmann equation 

The deviation function g has to be found by solving a kinetic equation. For weak enough 
interactions, quasiparticles remain well defined, and the problem reduces to solving a 
semiclassical Boltzmann equation for /^^^ in the presence of a static driving field: 



1 



/is|V5| 



Afe 



dt 



drive 



cr 
— ( 
T 



1 



/Xs|V5| 



crAfc 



dt 



coll 



cr'/' 



(2vr) 



C^,Mk,kr9^n'P^[Cglrk- (12) 



Here we have defined the driving term -D^^g, with which the Boltzmann equation can 
be written compactly as D = Cg. The kernel of the collision integral, C(jiyi>{k, k') 
defines the linear integral operator C, which has been described in detail in previous 
works concerned with thermal and electrical transport [TOl 120] . The form of C can 
be easily obtained from those references, taking into account that in the present case 
the dependence of ^ on cr requires to retain the cr, cr'-dependence of C. The two main 
ingredients of the kernel C are the occupation factors of in- and out scattering particles 
and the transition matrix elements, which enter the Fermi Golden rule for transition 
rates. At the level of the Born approximation, i.e., at the lowest order in the coupling 
a, the relevant transitions are depicted in Fig. |2j 

As we assume weak interactions, we can neglect the potential long ranged 
contributions to the collision integral due to the exchange of collective modes [27] . The 
latter are known to contribute to thermal transport, but they are expected to affect spin 
and electrical conductivity at best at subleading order in the interactions. 



2. 3. Variational principle for the spin conductivity 



From Eqs. ( [ll|l2| ) one sees that the solution of the Boltzmann equation satisfies 
as = 2T{g,D). This allows one to derive a convenient variational characterization 
of the spin conductivity (cf. [281 129] ) 



ex.. 



2T/i^ max 

g 



{~9, C~g) 



with {g, h) = ^ j 



d^k 



crlk' 



(13) 



Using the positivity and Hermiticity of the collision kernel. 



one can show that the 



solution of the Boltzmann equation maximizes the ratio in ( 13 ) 



The only dimensionless parameters of the problem at hand are a and the ratio /i/T. 
Thus, apart from the quantum unit of /i^/^, the spin conductivity is solely a function of 
these two parameters, cr^ = nl/hip^fi/T; a). In the remainder of this paper we compute 
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Figure 2. (a) The four classes of lowest order tree level Coulomb scattering processes, 
which enter the collision integral at leading order in a. The labels A, B stand for 
quantum numbers (ct, a, A), A being the "antiparticle" of A with all indices opposite. 
Apart from a factor which depends on the external momenta fci....4, the matrix elements 
are proportional to Ia^b, but independent of the indices a, a. Under the exchange of 
indices (a, I) — >■ {I, a) the two interfering diagrams (iii) and (iv) are multiplied by the 
same factor ±1, so that the square of their sum is invariant at the level of the Born 
approximation. At the Dirac point the occupation factors = /? are also invariant. 
Hence, the collision operator acquires an accidental SU{8) symmetry (if valley degrees 
of freedom are included) to leading order in the coupling a. (b) This does not hold 
anymore to order 0{a^) in the scattering rate. This is illustrated by the interfering 
diagrams (v), whose sum squared is not invariant under {a, I) — )• {I, a) if IaIb ^a^b- 

the scaling function ip, its minimum at //./T = and its asymptotic behavior at high 
chemical potential, ji/T ^ 1. 

3. Solution of the Boltzmann equation at the Dirac point 

Graphene has four species of massless Dirac quasiparticles, (z^, o") = (±1,±1), each one 
coming as particle (conduction band) and hole (valence band) excitations (Z = ±1). 
At the Dirac point, these 8 types of excitations are all equally populated by a finite 
temperature. The corresponding SU (8) symmetry of the non-interacting Hamiltonian 
is reduced by the Coulomb interactions to a SU{A) x U{1) symmetry. However, as we 
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will see below, at weak coupling this reduction of symmetry only manifests itself beyond 
the leading order interaction effects in the transport properties. 

Instead of working with valence band electrons with quantum numbers (cr, A, k, = 
\hvF\k\) and deviation function g^^j;, it is convenient to introduce hole excitations 

of positive energy with quantum numbers (5" = Xa,\,k = Xk,ek = Xek) and group 
velocity v = deu/dk = dlk/dk = k/k. [ji] Further we introduce the deviation functions 
aik ~ ^^(Tik' ^"^^ h^i{k) = lg„i{k). Positive /i,^^ corresponds to an excess of positive 
energy excitations as compared to equilibrium. A corresponding conjugation of the 
collision operator and the driving term lead to the Boltzmann equation in the form 
D = Ch, where 

C^i;a'i'{k,k') = lC^,y,{k,k')r 

b --. = ID (14) 

In the case of spin and electrical conductivity (where we define e as E/\E\) the driving 
terms are 

^:,s=*^/^(i-^°)- -f^!,. = - /?). (15) 

respectively. Note that at the Dirac point, in the absence of a background magnetic 
field, 

ft = i (16) 

' eMhvFk/T) + l 

only depends on k, since fi = B = 0. Thus the driving terms transform into each 

alk 

other under the exchange of indices {a, I) — )■ (/, a). Formally, 

& = TD^, (17) 

where we define the operator T with T^-i yi/ = S^^iiSiy. 

For electrical transport, it has been shown in Ref. ([261 HD]) that at the Dirac point 
the solution of the Boltzmann equation takes the form 

gliik) = g'^'ik) -> /.?,(fc) = A(7^^(fc). (18) 

From the symmetries of the Hamiltonian and the driving term, it follows similarly 
that in the case of magnetic driving the solution must take the form 

g:i{k) = ag'''{k) ^ ht.ik) = ^g'^'Ck)- (19) 

In both cases the deviation from equilibrium, h'^''^ only depends on the sign of the 
coupling of the quasiparticles to the external driving field, but is otherwise independent 
of the quantum numbers A, a, and u. 

A detailed calculation of the spin conductivity shows that within the Born 
approximation the two solutions at the Dirac point are actually identical, g^'^ik) = 

II The valley index v changes under such a transformation, too, i) — Xv. 
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g^'^{k). As a consequence, the dimensionless charge and spin conductivities are the 
same within the Born approximation (i.e., up to relative corrections of order 0(q;)), 

This result may come as a surprise since the Coulomb interaction, not being invariant 
under T, breaks the SU{8) symmetry of the non-interacting Hamiltonian and hence 
should lead to different values of the transport coefficients. Indeed, we will show below 
that and a* differ at higher orders in the coupling strength a. 

The reason for the accidental equality at the level of the Born approximation can 
be understood as follows. Upon exchanging simultaneously all quantum numbers A and 
a of incoming and outgoing particles, the transition amplitudes entering the collision 
integral change at most their sign. Applying Fermi's golden rule at the lowest order of 
perturbation theory in the collision integral, one finds that only matrix elements with 
identical transformation behavior interfere with each other, cf. Fig. [2j The common 
sign thus disappears upon squaring, and the resulting collision operator turns out to be 
invariant under the exchange of A and a, 

rp^CBov-n^rp—l ^(Born) (21) 

With Eq. Anh this immediately entails the relation h'^ = Th^ between the solutions, and 



thus the claimed equality of g^'^ and g^'^, and in particular Eq. (20). 

It is quite clear that to higher order in a transition matrix elements with different 
transformation behavior under T do interfere, which lifts the invariance of the collision 
operator found to lowest order. This is illustrated in Fig. [2] (b), which shows two 
interfering transition amplitudes whose contribution to the collision rate is sensitive to 
the exchange of / and a at order O(a^). 

A quantitative calculation of such higher order effects is very involved due to the 
large number of diagrams to be included in the collision operator. This is beyond our 
scope here. However, it may be interesting (and much simpler) to calculate the difference 
of spin and charge conductivities in examples of strongly coupled critical theories, which 
can be exactly solved by the AdS-CFT correspondence [3U1 ED 1321 E3] • An experimental 
measurement of both conductivities and in particular their difference, might provide an 
interesting test for the presence of strong interactions. On one hand such an experiment 
could quantitatively compare the two conductivities, which are both of the order of their 
natural quantum units. On the other hand it could potentially disentangle interaction 
effects from disorder effects, which so far dominate the existing experiments 

4. Spin conductivity as a function of chemical potential 



In this section we work with units in which = h = vp = 1, reinserting them only in 
final results. 
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4-1- Separation of relaxation time scales in two dimensions 

Scattering problems in two dimensions have the interesting property that small angle 
forward scattering is strongly enhanced, without destroying immediately the Fermi 
liquid, unlike in one dimension [311 CH]- This phenomenon introduces a separation of 
time scales between the relaxation which establishes equilibrium among quasiparticles 
with equal direction of motion, and the relaxation to global equilibrium. This feature 
was exploited in previous work, which mostly focussed on the critical window /i < T, 
see Refs. [TOl EHl [35]. Since, the respective timescales are parametrically different at 
weak coupling, one can obtain an approximate, but asymptotically (in small a) exact 
solution of the Boltzman equation by projecting it onto the slow relaxation modes. 
At the Dirac point this approximation is parametrically controlled by the logarithm 
of the weak coupling strength. As we will show below, the approximation becomes 
even much better at finite chemical potential for weakly screened Coulomb interactions, 
which strongly enhances small angle scattering and the tendency to establish partial 
equilibrium among particles with equal group velocity. 

The scale of fast relaxation rates of generic deviation profiles gai{k) can be estimated 

as 

where the denominator {g\fj^f,{l - f^^M = J2ai I 4^9^ikffk{'^ - fik)9aik normalizes the 
collision rate. The collisions at /i < T are dominated by forward scattering events. 
Their logarithmically divergent cross section is only cut-off by life-time broadening of 
the quasiparticles or higher order corrections in the interactions. (TU] This leads to a 
scattering rate 

^r^lil^ <T)- {V,TfT\og{[V,Tf) (23) 

where Vq = V{q ~ T) is a typical interaction strength at thermal momentum 
transfers [I^ (VqT ~ a for Coulomb interactions, while for short range interactions 
^(9) = ^! independently of q). In the gated case, 3> T, the fast relaxation is still 
dominated by forward scattering. Taking into account kinematic constraints, it can be 
estimated as an integral 

hr-^,{^, » T) ~ - r dct>—^ {[^^V{q = 2/isin(0/2)]^ + [^^V,f) ,(24) 

The first term refers to scattering processes as illustrated in Fig. [3| cj) being the scattering 
angle between incoming and outgoing momenta, ki and ^4. The second term is due 
to small angle scatterings which essentially preserve momenta, now being the angle 
between incoming momenta ki^2- The factor l/|sin(0)| is crucial and originates from 
two-dimensional kinematics and phase space constraints. The logarithmic divergence 
at = 0, when all particles move parallel to each other, is only cut off by lifetime or 
screening effects, while the divergence at = tt is regularized at angles |vr — 0| ~ 
For short range interactions the integral contributes a factor (/iVo)^ and an additional 
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logarithm \og[n/ 5E), where 5E is the hfe time broadening of quasiparticles. However, 
for weak Coulomb interactions and T / ^ < a, the integral is strongly dominated by small 
momentum transfers g ~ go ^ A^, leading to 

^^f;st,cb(-" » T) ~ — log(«/i/5E). (25) 
/i 

The strong forward scattering leads to a rapid equilibration of quasiparticles moving 
in the same direction. The remaining slow modes correspond to deviation functions g, 
where each subset of quasiparticles with the same is already in equilibrium. The only 
such function g^i{k) which is also compatible with linear response (azimuthal dependence 
g e- V and spin dependence ^ ~ cr) is 

9lik = ^ ^^'^ik- (26) 

Indeed, the quasi-equilibrium along fixed angles suppresses small angle scattering, the 
slow relaxation for the gated case being estimated as 

K^K-" » T) ~ - rd<p'^^^^ [f,V{q = 2/isin(0/2)]^ , (27) 
P Jo I sini^</)j| 

where the extra factor sin^(0/2) arises from the fact that quasiparticles traveling at 
similar angles are nearly equilibrated (see also Eqs. ( 30p2 ) below). For short range 



interactions this rate scales as (T^//i)(Vo/i)^, which is logarithmically smaller than the 
fast relaxation rate. For weak Coulomb interactions the ratio between the slow rate 



^^swcb(/^ » T)r^a'— log(l/a), (28) 



and the fast rate (25) is even smaller, of order O(a^), apart from logarithmic factors. 
This slow relaxation is dominated by head-on scatterings between opposite spins with 
small momentum transfer, see Fig. [sj Likewise, close to the Dirac point, hT~^^^{iJ, < 



T) ~ (VoT)^T is smaller than the fast mode (23) by the logarithm which is absent in 
the slow relaxation. 

4-2. Spin conductivity from the slow mode approximation 

The ratio between slow and fast relaxation rates is a small parameter which allows us 



to obtain an excellent estimate of the spin conductivity by injecting the slow mode (26) 



into the maximum principle (13). Equivalently, we observe that the separation of time 
scales assures that the solution of the Boltzmann equation is proportional to the slow 
mode up to parametrically small corrections. 

This principle has been used in Ref. [TD] to obtain the electrical conductivity at 
fi = within a weak coupling approximation up to logarithmic corrections. As shown 
above the same numerical value applies to the spin conductivity, and thus 

^. - (29) 

The separation of time scales also allows for a parametrically controlled derivation of 
hydrodynamic equations from the semiclassical Boltzmann equation. 
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Figure 3. The dominant scattering channels in the degenerate Fermi hquid regime 
S> T) consist in (a) head-on colHsions with fci + ^2 ~ ^3 — ^4 ~ 0, and (b) scattering 
events with nearly identical pairs of incoming and outgoing momenta (fci, ^2) ~ (^3) ^4) 
or (kijk^) [not shown]. For generic distribution functions g small angle scattering 
q II dominates the scattering cross-section. The relaxation of the slow mode, 
g oc ak/k, proceeds essentially via the two shown channels, involving scattering 
between opposite spins cti = —<J2- 



4-3. Degenerate Fermi liquid regime ^ T and 2d metals 



To leading order in /i/T the collisions are dominated by scattering of opposite spins, 
(^i,t)^2,i) — ^ (^45!) ^3) 4); either with ~ ki,k4 ~ /c2 or with ki ~ —^2,^3 ~ — ^4 
(head-on collisions), see Fig. [s] These are the only momentum and energy conserving 
processes which are not strongly suppressed when the deviation function has relaxed to 



the slow mode, g^^f: — )■ g^^. In that case the collision matrix element in (13) simplifies 



to 



xF°{ki,k2,k3,k4 



dki dk2 dfca 
(27r)2 (27r)2 (27r)2 

T++-\-+{ki, k2, k^, k^ 



S{\ki\ + I/C2I - Iksl - \k4\) 



where F° is the product of occupation factors. 



F'>{ki,k2,k3,k4) 



1 



1 



1 



k2 ks 



1 



(30) 



(31) 



+ le^' ^/^Igv^ --' + leA^ ■^^' + 1 
while the matrix element T+_|__(_+ was defined in Eq. ([s]). The last term in (30) arises 
from the "transport factor" [^"(^1,0") -|- g^{k2,—a) — g^iji^.—a) — ^°(/c4, cr)]^. In the 
limit fi ^ T, the integral (30) can be evaluated analytically by observing that the last 
two factors in the integrand only depend on the scattering angle between ki and k^, 
while the integrals over the other variables can be carried out analytically. The two 
dominant scattering channels shown in Fig. [3] both contribute an equal amount to the 
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matrix element, with the result 



7.0\ 



37r2 



2n 



COS^ 0/2 . 2 



sm 



sm( 



fiV q = 2fism 



.(32) 



For weakly screened Coulomb interactions, cf. Eq. ([T]), this yields for T / ji <^ 



a 



^3 



3 

2a'T^ 



p2TT 


, cos^ 0/2 


sin^ 0/2 




^ \ ■ A\ 
sm0 


(2sin0/2 + 4a)2 


j.n( 


^2^) 


+ 0{a) 





(33) 



The matrix element of the slow mode with the driving term evaluates simply to 
T 



{g^^D) = - In [2 cosh(^/2r)] 



7[ ' 27r 
We thus find the maximum of the spin conductivity functional as 

.cb,,... ^,_o^ Cg'^Dy _^.l{^i/Tf ^ 



a'^^'i^i > T) = 2T^ 







-1 















(34) 



(35) 



up to small relative corrections of order 0{T/fi). The corresponding result for short 
range interactions tends to the /i-independent limit. 



(/^ » T) 



(36) 



We confirmed numerically, by extending the maximization of the variational 
principle to a larger function space g, that g^ indeed gives an excellent approximation 
to the actual maximum of (13). 

For small and moderate values of fi/T, {g^, Cg^) has to be be evaluated numerically. 
The resulting full scaling function for the spin conductivity is plotted in Fig. |4j We also 
display the spin diffusion coefficient Ds, which is related to as by the Einstein relation 
Ds = cTs/x where x = ~ (fe»J)2 lii[2 cosh(/x/2T)] is the Pauli susceptibility. At the Dirac 
point Dg ~ 1/T, while in the degenerate Fermi liquid regime one finds ~ fJ'/T'^- 

2 

The scaling as ~ ^(/i/T)^ is easy to understand with similar arguments as those 
given in the introduction for the conductivity at the Dirac point. The scaling is expected 
from a Drude-type estimate as ~ /i^nrinei/?^eff, where mes = Pf/vf = 1^1 v\ is the 
effective inertial mass and n ~ {fi/hvpY the density of carriers, while Tinei ~ hfi/T'^ 
is the inelastic transport scattering time of a 2d Fermi liquid. We stress that the 
transport time coincides with the slow relaxation time, which in two dimensions is 
parametrically smaller than the collision time, which scales like Tfast- The latter controls 
the quasiparticle life time, but not the transport. In contrast to the collision and life 
times, the transport time does not contain an extra logarithm of fi/T, a fact which was 
missed in earlier studies of quasi-two dimensional metals |36j . 

In the case of screened Coulomb interactions the transport time is enhanced by 
a factor log(a) which reflects the enhanced scattering at small momentum transfer. 
Apart from this logarithm the transport time is independent of the bare strength of the 
interactions a. 
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Figure 4. Spin conductivity as and spin diffusion constant TDg as a function of 
chemical potential. The diffusion constant behaves as Dg ^ 1/T close to the Dirac 
point, and scales as ^ |/i|/T^ in the degenerate Fermi regime |^| ^ T. The inset on 
the left shows the approach of a-g to the analytical prediction Eq. ( 35 ) (dashed line). To 
regularize small momentum transfers Thomas-Fermi screening has been incorporated 
assuming a coupling a = 0.01. Note that the screening slightly increases the value of 
asifJ. = 0). 



5. Discussion and Conclusion 

The minimal spin conductivity of the order of the quantum unit fi^/h reflects the 
strong inelastic scattering among quasiparticles at finite temperatures close to the 
Dirac point. It is one of several interesting hallmarks of the strongly coupled electron- 
hole plasma, together with a minimal collision-dominated electrical conductivity and 
a minimal viscosity to entropy ratio at /i = 0. All these aspects are very similar to 
the phenomenology near quantum critical points, and arise here due to the marginal 
irrelevance of the Coulomb coupling. However, the disorder free, collision-dominated 
spin conductivity is also well-defined away from the Dirac point, as well as in any clean 
metal without spin orbit coupling. The analytical results we obtained for graphene in 
the degenerate Fermi liquid regime are general and immediately carry over with minimal 
adjustments to the case of 2d metals with simple Fermi surfaces. 

The simplest way to observe collision-limited spin diffusion consists presumably in 
preparing a local spin polarization by an external field, and monitoring its spreading 
(with diffusion constant Dg) after switching off the field. In reality this diffusion will be 
limited in space and time by a finite (even though rather long) spin flip time. 

Alternatively, the following Gedankenexperiment involves a spin current source to 
measure ag in the analog of an electrical two-point measurement. Imagine to contact 
the sample by the two opposite edges of a spin Hall conductor [37] , in which a spin Hall 
voltage is induced. This spin-Hall source then injects oppositely polarized currents into 
the sample, but no net current. The resulting total spin current in the sample will be 
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controlled by the collision-limited spin conductance of the sample (proportional to as, 
the width and the inverse length of the sample) and the contact resistances. 

We hope that future experiments in suspended graphene with unscreened Coulomb 
interactions will explore these interesting facets of collision-dominated transport. 
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